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Abstract 



We have proposed a quantum model for the magnetic multi- valued recording in this pa- 
per. The hysteresis loops of the two-dimensional systems with randomly distributed mag- 
netic atoms have been studied by the quantum theory developed previously. The method 
has been proved to be exact in this case. We find that the single-ion anisotropies and the 
densities of the magnetic atoms are mainly responsible for the hysterisis loops. Only if the 
magnetic atoms contained by the systems are of different (not uniform) anistropies and 
their density is low, there may be more sharp steps in the hysteresis loops. Such materials 
can be used as the recording media for the so-called magnetic multi- valued recording. Our 
result explained the experimental results qualitativly. 
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Magnetic thin films with perpendicular "easy-axis" anisotropy have attracted much 
attention these years both experimentally and theoretically since such systems have good 
potential to be taken as recording media In practice, many efforts have been 

made to increase the recording density of the devices, and several proposals have been 
made to achieve this end. However, the recording density has already come to a limit of 
the conventional scheme so that one must find new approaches. Recently, experimental 
studies have been made on the magnetic multi-valued (MMV) recording which is believed 
to be the next strategy of high density recording The key point is that there 

should exist more metastable states which are stable enough to record a message. Thus 
the media for MMV recording must possess several sharp sub-steps in its hysteresis loop. 
Experimentally, such a phenomenon had been confirmed in some kind of magnetic layered 
systems [jlj -B. However, the theoretical origin is not yet very clear. 

On theoretically side, a quantum theory for the coercive force of a magnetic system || 
has been developed based on some previous works In such a quantum approach, 

the concept of the metastable state was adopted, and the magnetic excitation gap was 
defined to be the order parameter to monitor the stability of such metastable states. The 
coercive force can be determined by the condition that the gap comes to zero ||. The 
so-called "capping effect" in a double-film structure || has been explained successfully 
by the quantum theory ||. 

The present work is devoted to proposing a theorectical explanation for the MMV 
recording which is confirmed by a randomly magnetic thin film 0. We first propose 
a model Hamiltonian for such a system in which the magnetic atoms are distributed 
randomly , then studied its hysteresis loop by the quantum method which is proved 
to be exact in such case. The results are averaged for samples finally to overcome the 
fluctuations of the distribution. We show that: 1) When the magnetic atoms contained 
by the system are of uniform single-ion anisotropy, there should be only one sharp step 
in the hysteresis loop. 2) When the magnetic atoms are of different anisotropics but the 
density is higher than a critical value, the hysteresis loop is highly smoothed and there 
are no obvious steps. 3) only when the magnetic atoms are of different anisotropies and 
their density is low, more sharp steps can be clearly observed in the hysteresis loop. 

The Hamiltonian is given as: 

H = \ll • Si - E A(Sf) 2 -h^s; (i) 

i,j i % 

where { % } are the lattice sites randomly distributed in the x-y plane in which the 
magnetic atoms are occupied. {Di > 0} are the single-ion anisotropies. An external field 
is applied along the z axis. Jij = J is the exchange constant and only the nearest-neighbor 
interaction is considered. In real magnetic materials, the single-ion anisotropic constant 
D is usually much smaller than the exchange constant J. 

Following Refs. H,||, a local coordinates (LC) system (xi,yi,Zi) can be introduced 
to the Hamiltonian. The Xi and Zi axes in the LC system are rotated by an angle 0j 
which may be different from site to site, while the jji axis is not rotated. It is helpful to 
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apply a Bose transformation such as Holstein-Primakoff (H-P) H], Dyson-Maleev (D-M) 
]l(J-[]lT] and the complete Bose transformation (CBT) [12| to study the spin systems. In 
a harmonic approximation, these transformations are the same. However, we will use the 
CBT in this paper because it can present high-order terms correctly. 
After the LC transformation and the CBT, the Hamiltonian becomes: 



where 



H = U + H 1 + H 2 + ■ ■ ■ (2) 



Uo = const. — 2 X] cos (^j ~~ ~~ h 2J cos 9iS% 

+ ^A(25 i -l)5 i cos 2 i (3) 
H 2 = E^(%^ + E G mWK + <+^')- (4) 



The coeffiecients are 



F u (9) = ^ (sin 2 0,-2 cos 2 0,) 



+ SjJij cos(9i — 9j) + h cos 9i, (5) 
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FiM = -^^V^-tl + cos^-^)], (6) 



= -^2^(2^-1) A sin 2 f , (7) 



= -^V^-tl-cos^-^)], i^j. (8) 

{0i} can be obtained by minimizing the ground state energy: dUo/d9i = 0, which yield: 
JijSj sin(9i - 9j) + h sin 9 % + A (25* - 1) sin 0, cos ^ = 0. (9) 

i 

Equations above are the same as the condition of Hi = 0. H 2 can be diagonalized by a 
Bogolyubove transformation. So, we have: 

H~U^ + Y f e i ata i + --- (10) 

In general, the nonlinear equations (||) m ay have many solutions in a defmate external field 
h. For each solution, the spin state may be a metastable one only when the excitation 
energy {e^} based on such a solution are all positive. In another word, the excitation 
should have a positive gap A(h) > 0, otherwise, such a spin state is not stable. So, 
if every metastable states of the systems have been investigated when one altering the 
external field from positive to negative, the hysteresis loop can be obtained while the 
system transits from one metastable states to another one. 



3 



Before carrying out a general investigation, we will first study some particular systems 
which are illustrated in figure 1 in order to get some informations. In figure 1 and 
other figures, we suppose that the magnetic atoms represented by the circles have the 
anistropy D\j J = Di(2Si — 1) /J = 0.1, and those by the triangles have the anistropy 
D 2 / J = Di(2Si — 1)/ J = 0.02. It is also supposed that Si = S 2 for convenince. It should 
be noted that such systems have a character that the magnetic atoms are all coupled 
together since we have supposed that only the nearest-neighour interaction is considered. 

Eqs. may have two kinds of solutions: (1) The trivial solutions: {#j = 0, or 7r}. 
(2) The non-trivial soltions: {9i ^ 0, or tt}. The number of solutions is certainly very big, 
and any of them may be a metastable state if the elementary excitations have a positive 
gap. However, after the numerical calculations for the systems which are illustratated 
in figure 1, we find that only two solutions among so many ones can be metastable 
within a definite region of external field. They are: (1) {9,i = 0} (2) {6i = 7r}. Other 
spin configurations are all unstable. Of cause, we can not give a vigorous proof for that 
in a general case, but one can understand the above numerical results as follows. First, 
the "easy-axis" of each magnetic atom is along the z axis so that the spins all prefer 
to paralell or antiparalell the z axis, and there are no such atoms which have "in-plane 
easy-axis" or "easy-plane" anisotropies, as the result, the non-trivial spin configuration 
{9i ^ 0, 7r, i — 1, 2, • • •} is not likely to appear. Second, since the exchange interaction J is 
usually much larger than the single-ion anistropy D in real material, the spins are willing 
to paralell with each other. Since the spins are all coupled together in the systems studied 
(figure 1), such spin configurations that some spins are up while others are down can not 
be stable. This is an explanation rather than a proof, but at least, one may understand 
the numerical results that only two metastable states may exist for a system in which the 
spins are coulped together. 

Applying the solution {8i = 0} into the Hamiltonian, we find that the excitation energy 
{ej} can be obtained by diagonalizing the matrix {Fij} since the matrix {Gij} is zero. In 
the Appendix, we prove that {e/} are the rigorous excitation energies of the Hamiltonian 
in this case although they seem to be obatained by a harmonic approximation. Thus we 
get exactly 

A(h) =Min[e i ], (11) 
and the coercive force of the model can be determined by A(h c ) = 0. 



Following are the numerical results we obtained. 

Case 1. Only kind of magnetic atoms are distributed in the systems: Si = S and 
(2Si — l)Di = D. In this case, the matrix {F^} can be rewritten into the following form: 

{F lJ }\ di=0 = D + h + {F: j } (12) 

One can check at once that the matrix {F(j} has the lowest eigenvalue 0. Then the gap 
will be: 

A(h) = D + h (13) 
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The coercive force of the system can derived: h c = —D. 

Thus, no matter how the atoms are distributed in the pattern, if there are only one 
kind of magnetic atoms, the coercive forces are always the same. In practice, one makes 
the systems amorphous to minimize the domain's size to improve the recording density. 
The above results can warrant that every domains possess the same coercive force so that 
they can be used in a same way. 

Case 2: The system contains more than one kind of magnetic atoms. The results for 
the systems shown in figure 1 are listed in Table 1, from which we can find that the 
coercive forces h c are strongly dependent on the distributions of the magnetic atoms. 
Generally, the coercive force turns smaller when the samll-anisotropy magnetic atoms are 
relatively increased in the system. But the results are different for the systems in which 
the magnetic atoms are coupled differently although the numbers of the two kinds of 
atoms are the same. 



We will discuss the general behavious of our random model's hysteresis loop. Through- 
out this paper, a 10 x 10 lattice is discussed just to illuminate the main physical idea, and 
the samples are avaraged in the end to take account of the fluctuations. 

A typical distribution is illustrated in figure 2. One may find that there are many 
so-called "Isolated Islands" (II) in the patterns. Within each II, the magnetic atoms 
are coupled toghter so that they can be treated following the method mentioned above. 
Howevr, the couplying between those lis is zero. So, such lis must be able to be considered 
independently. Actually, the matrix F^ has the following form: 

/ [ C« ] 

[ ] 

[ ] 

K [ cw ] 

where [C^] are the sub-matrixes for each lis. They are found not to couple with each 
other in the matrix. Thus, the eigenvalues of the matrix {Fij} must be those of the 
sub-matrixes [C^]. Without loosing any genarality, it can be supposed that the coercive 
forces h c of all the lis have been arranged from small to large: \h % c \ < \h 3 c \,i < j. 

The hysteresis loop of the entire structure can be obtained as follows. As the external 
field h reaches hi, the minimum eigenvalue of the sub-matrix [C^] is zero while that in 
others are still positive, that means the first II is not stable while others are still stable. 
It has been discussed that each II only has two metastable configurations: the spin-up 
and the spin-down configuration. So, the first II will turn to its spin-down configuration 
in the vicinity of hi while the other lis will remain in there spin-up configuration. Fol- 
lowing the same reason, the ith II will turn to its spin-down configuration while other 
IPs configuration are not changed as the external field h reaches the ith coercive force h l c . 
If we do that step by step, we can obtain the hysteresis loop finally. For instance, the 



(14) 
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hysteresis loop for the system illustrated in figure 2 has been shown in figure 3, where 
there are sharp decrease in the vicinites of the coercive forces h l c . 

However, we must do the sample avarages to overcome the fluctuations. In our calcu- 
lations, 1000 samples are averaged for a definite case. The results are presented in figures 
4-5. In the case that the density of magnetic atoms is 40% and the two kinds of magnetic 
atoms are equally sputtered, we average two groups of 1000 smaples to get two final re- 
sults. They are compared in figure 4. However, one may find that the difference of the 
two lines even can not be detected which means the number of the samples for average is 
large enough to overcome the fluctuation. Figures 5a-5d present the hysteresis loops for 
the systems in different cases. From figures 4, 5a-5c, one may find that more sharp steps 
can be apperantly found in the hysteresis loop. Futhermore, when the density is lowered, 
the steps are sharper and more sub-steps in the hysteresis loops may appear (figure 5a). 
This can be understood as follows. When the density is low, the possiblities of apperaing 
some definite structures will be high. As a result, there may be a distint decrease of the 
magnetization in the vicinity of the coercive force for such a structure. Actually, a limit 
case is that there are only two differnt atoms in the lattice, in this case, the possibility 
for the distribution that the two atoms are seperated should be greatly larger than that 
they are coupled together. Thus, they may be a very clear multi-step shaped hystersis 
loop. 

When the density of the magnetic atoms is high, especially when the density is larger 
than the percolation value, the hysteresis loop is greatly smoothed and the step is almost 
undetectable (figure 5d). Actually, near the percolation value, the distributions of the lis 
are quite complicated, and any pattern is possible. Since the coercive force of each II is 
strongly dependent on the distribution, any value of the coercive force is then possible 
to appear. Thus, there should be an infinite number of metastable states which are all 
different in the system. This is very similar to that in the spin-glass system, although 
such a system is quite different with that one. So, one should not use such materials for 
recording. 

To summarize, in this letter, we have investigated the hysteresis loops of the 2-d 
systems with randomly sputtered magnetic atoms. The method is proved to be rigorous 
in this case. The results show: a multi-step shaped hystersis loop can be achieved only if 
the following two conditions can be satisfied: 1) the sputtered magnetic atoms must have 
different coercive forces. 2) the density of the magnetic atoms must be lower than the 
percolation value. Such materials have the potential to be considered as the recording 
media for MMV recording. If the first condition is dissatisfied, there is only one sharp 
step in the hystersis loop; if the second condition, there are no sharp step in the hystesis 
loop and such material can not be used for recording. 
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Appendix 



In this appendix, we will prove that {e^} are the exact excitation energies of the Hamil- 
tonian (1) in the case of 0j = 0. 

Since &i = 0, the Bose transformation can be applied naively to Hamiltonian (1). 



Following Ref. ||12||, the CBT is given as: 



1=0 

oo 



1=0 

oo 

1=0 



l a t l(l i +1 



(15) 



{A}}, {Bj} are the coefficients of Bose expansion which are dependent on S. We can also 
derive the expansions for the single-ion anisotropy terms: 



{SI? = Y,G\afa l 



(16) 



1=0 



Applying the CBT to the Hamiltonian, we have a transformed Hamiltonian H which 
has exactly the same eigenvalues as the Hamiltonian (1): 



H = U + H 2 + H 4 + 



where 



Ho 



(17) 

(18) 
(19) 



Of course, Hamiltonian H is still impossible to solve exactly However, some eigem- 
states can be obtained exactly. H2 can be diagonalized by a orthogonal transformation: 



(20) 



After the transformation, we find 



H2 — E € n a n a n 



(21) 
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where {e„} are the eigenvalues of the matrix {F^} and P mn can be found from calculating 
the eigenvectors of the matrix {-F^}. 

From the CBT, we find that every terms in the remainder interaction Hj = H4 + Hq + 
■ ■ ■ contains equivalent numbers of the creation operators a + and annihilation operators 
a and contains at least two annihilation operators. For example, H4 is found to be 

H A = J ii l A o A 2 a t 2a2 i + A o A i dj 2 a 2 + M. A i a t a f a i a j 

+ C\Ciata+{al + a))] -J]A« - hJ2 A WX- (22) 

i i 

Then, it is easy to prove that 

H I ~a+\0) = J2Pn, m Hia+ l \0) = (23) 

m 

As the result, 

ff(5+|0»= CB (5+|0» (24) 

So, such one-magnon eigenstates a n \0) are the exact eigenstates of the system, and {e^} 
must be the exact excitation energies of the system. 
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Captions: 

Figure 1: Some particulr pattterns studied in this paper 

Figure 2: A typical distribution of the magnetic ions in the system we sudied 

Figure 3: The hystersis loop of the system illustrated in figure 2. 

Figure 4: Comparison of the hysteresis loops obtained by averaging two 
groups of 1000 samples seperately 

Figure 5: Hysteresis loops for the systems in different cases 
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